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Abstract 

Intermediate energy heavy ion collisions open the unique possibility to explore the 
Equation of State {EOS) of nuclear matter far from saturation, in particular the 
density dependence of the symmetry energy. Within a relativistic transport model 
it is shown that the isovector-scalar 5-meson, which affects the high density behav- 
ior of the symmetry energy density, influences the dynamics of heavy ion collisions 
in terms of isospin collective flows. The effect is largely enhanced by a relativistic 
mechanism related to the covariant nature of the fields contributing to the isovector 
channel. Results for reactions induced by ^'^'^Sn radioactive beams are presented. 
The elliptic flows of nucleons and light isobars appear to be quite sensitive to micro- 
scopic structure of the symmetry term, in particular for particles with large trans- 
verse momenta, since they represent an earlier emission from a compressed source. 
Thus future, more exclusive, experiments with relativistic radioactive beams should 
be able to set stringent constraints on the density dependence of the symmetry 
energy far from ground state nuclear matter. 
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The high density behaviour of nuclear symmetry energy Egym^PB) is very 
important for understanding many interesting astrophysical phenomena, but 
it is absolutely not constrained by the predictions from several relativistic 
and non-relativistic [1,2] models of nuclear matter. The results can be roughly 
classified into two groups, i.e., one where the EgymiPB) rises, Asy-stiff, and one 
in which it falls with increasing density, Asy-soft [3]. An increasing EgymipB) 
leads to a more proton-rich neutron star whereas a decreasing one would make 
it more pure in neutron content. As a consequence the chemical composition 
and cooling mechanism of protoneutron stars [4,5], mass-radius correlations 
[8,9], critical densities for kaon condensation in dense stellar matter [6,7] as 
well as the possibility of a mixed quark-hadron phase [10] in neutrons stars 
will all be rather different. It has recently been argued by means of simple 
thermodynamics considerations that the onset of the quark phase has a strong 
sensitivity to the behaviour of Egym^pB) even for not very large asymmetries 
[11]. 

The search for Esym{pB) around saturation density has driven a lot of theo- 
retical and experimental efforts. It seems to be rather well established that 
heavy-ion collisions {HIC) at cyclotron energies can give the possibility to 
extract some information on the symmetry term of the nuclear Equation 
of State (EOS) in region below and/or sligthly above the normal density 
[3,12,13,14,15,16]. On the other hand it is quite desirable to get information 
on the symmetry energy at higher density, where furthermore we cannot have 
complementary investigations from nuclear structure like in the case of the 
low density behaviour. Indeed HIC provide a unique way to create asym- 
metric matter at high density in terrestrial laboratories. Calculations within 
transport theory show that HICs around lAGeV allow to reach a transient 
state of matter with more than twice the normal baryon density. Moreover, 
although the data are mostly of inclusive type (and the colliding nuclei not 
very neutron rich), quite clearly a dependence of some observables on charge 
asymmetry is emerging. 

In this paper we show that a relativistic description of the nuclear mean field 
can account for an enhancement of isospin effects during the dynamics of 
heavy-ion collision. In particular future experiments with radioactive should 
be able to provide information on the vector part of isovector mean fields from 
collective flow analyses. 

The isospin dependence of collective flows has been already discussed in a 
non-relativistic framework [12,17] using very different EOS with opposite be- 
haviours of the symmetry term at high densities, increasing repulsion {asy- 
stijf) vs. increasing attraction {asy-sofi). The main new result shown here, in 
a fully relativistic scheme, is the importance at higher energies of the micro- 
scopic covariant structure of the effective interaction in the isovector channel: 
effective forces with very similar symmetry terms can give rise to very different 
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flows in relativistic heavy ion collisions. 

We start from the Relativistic Mean Field (RMF) picture of the hadronic 
phase of nuclear matter [18] which has been extensively used to the study of the 
EOS for symmetric and asymmetric matter. The RMF describes an interact- 
ing system of nucleons (described as Dirac spinors) and meson classical fields. 
The most common treatment of the isospin dependent part of the interaction 
is based on the introduction of an effective p-meson field (vector-isovector) 
which can account for the known value of symmetry energy at normal density. 

However, a full description in a relativistic framework in principle should rely 
on the balance between a scalar and a vector field as stressed in some pa- 
pers within the Hartree approximation [20,19], and as naturally accounted 
for within the Dirac-Hartree-Fock {DHF) [21,22] or the Dirac-Brueckener- 
Hartree-Fock (DBHF) [23] schemes. One could argue that the scalar-isovector 
S meson (ao(980)) is too heavy or that the scalar field in the isoscalar channel 
may be due to the two pion correlation, while there is no equivalent possibility 
in the isovector channel. On the other hand we stress that the vector-isovector 
field in the RMF effective picture has not to be viewed as coming only from 
the exchange of the S meson. In fact in the nuclear system contributions to the 
isovector channel are mostly due to isoscalar mesons, as clearly shown by DHF 
or DBHF scheme [22,23,24,25], through important exchange and correlation 
terms. Therefore when in the following we will refer to 5 field indeed we mean 
a, 5 — like field, i.e. the scalar isovector part of the effective "interaction" . 

In recent years some efforts have been devoted to the effects of the scalar- 
isovector channel in finite nuclei, [26,27,28]. Such investigations have not shown 
a clear evidence for the 5-field and this can be understood considering that in 
finite nuclei one can test the interaction properties mainly below the normal 
density, where the effect of the 5— channel on symmetry energy and on the ef- 
fective masses is indeed small [19] and eventually could be absorbed into non 
linear terms of the p field. Moreover even studies of the asymmetric nuclear 
matter by means of the Fermi Liquid Theory [19] and a linear response analy- 
sis have concluded that some properties, like the borderline and the dynamical 
response inside the spinodal instability region, are not affected by the S field 
[29]. 

Here we show that heavy-ion collisions around lAGeV with radioactive beams 
can provide instead a unique opportunity to spot the presence of the scalar 
isovector channel. In fact, due to the large counterstreaming nuclear currents 
one may exploit the different Lorentz nature of a scalar and a vector field. 

All models, including RMF, allow to use the parabolic approximation for the 
description of the EOS of asymmetric nuclear matter: 

EipB,I)=E{pB) + E,,^{pB)P (1) 
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Fig. 1. Total (kinetic+potential) symmetry energy (in MeV) as a function of the 
baryon density. Solid: RMF- {p+5). Dashed: RMF-p. Short Dashed: RMF-Dp. 
In the insert the density behaviour of the p coupling, fp (in fm?), for the three 
models is shown. 

^^-^ is the asymmetry parameter. When both p-like and 5-like 
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channels are considered EgymiPB) can be written as [19]: 
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with E* = ykp + M*^, M* the effective Dirac mass and fp^s = {gp,s/'m'p,sy 
are the coupling constants of the isovector channels. 

We see that, when 6 is included, the observed 04 = Egymipo) value actually 
assigns the combination [fp — /^(^)^] of the {p,6) coupling constants, for 

futher details see ref.[19]. In Fig.l we report the density dependence of the 
symmetry energy for three different models: one including only the p field 
{RMF - p), the other with (p + S) fields {RMF - {p + 5)) and the last, 
{RMF — Dp) , including only p field but with a covariant density dependence 
of fp, see ref. [26,30]. This is tuned just to give at high density the same 
Esyra{PB) of the the {RMF — {p + 5)) case. As shown in the following, the 
latter is useful for disentagling in the reaction dynamics the effects due to 
the difference in Esym{pB) from those directly linked to the strenght of the p 
vector field. 

Thus these models parametrize the isovector mean field either by only the 
vector field with fp — 1.1/m^, or with a balance between a vector field with 
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f/j = and a scalar one with fs = 2.4/m^, or finally by a normal 

density coupling /p(po) = l-lfni^- but with an increasing density dependence 
as shown in Fig.l (insert). We stress again that in RMF — {p+d) the symmetry 
energy is coming from a balance between a scalar attraction, {S — like), and 
a vector repulsion, (p — like), which is now roughly three times larger than in 
the RMF - p case. 

The choice of fs is fixed relatively weU by DBHF [23] and DHF [22] cal- 
culations. Therefore the effect described in the following is not artificially 
enhanced, but based on a reliable estimate available at the moment. In any 
case the aim of this work is to present some qualitative new features expected 
in the reaction dynamics, in particular for collective flows, just due to the 
introduction of a scalar effective field in the isovector channel. 

Collective flows in heavy ion collisions give important information on the dy- 
namic response of excited nuclear matter [31,32]. In particular the proton- 
neutron differential flow F'^^\y) [17] has been found to be a very useful probe 
of the isovector part of the EOS sine it appears rather insensitive to the 
isoscalar potential and the in medium nuclear cross section. The definition of 
the FP^{y) is 

1 N{y) 

where N{y) is the total number of free nucleons at the rapidity y, p^,. is the 
transverse momentum of particle i in the reaction plane, and Tj is and -1 
for protons and neutrons, respectively. 

For the theoretical description of heavy ion collisions we solve the covari- 
ant transport equation of the Boltzmann type within the Relativistic Landau 
Vlasov (RLV) method [33] (for the Vlasov part) and applying a Monte-Carlo 
procedure for the collision term. RLV is a test particle method using co- 
variant Gaussians in phase space for the test particles. The collision term 
includes elastic and inelastic processes involving the production/absorption 
of the A(1232Mey) and A^*(1440Mey) resonances as well as their decays 
into one- and two-pion channels. Details about the used cross sections for 
all possible channels can be found in Ref. [34]. An explicit isospin-dependent 
Pauli blocking term for the fermions is employed. Asymmetry effects are suit- 
ably accounted for in a self-consistent way with respect for the RMF models 
discussed above. 

A typical result for the ^^'^Sn +^^'^ Sn reaction at 1.5AGeV (scmiccntral colli- 
sions) is shown in Fig. 2. The error bars are related to statistical fluctuations 
due to the Monte-Carlo nature of the simulations. 
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We notice that the differential flow in case of the RMF— (full circles and 

solid line) presents a stifTer behaviour relative to the /?Af F — p (open circles) 
model, as expected from the more repulsive symmetry energy Egym^PB) at 
high baryon densities, see Fig.l. On the other hand it is quite surprising that 
a relatively small difference at 2po can result in a such different collective flows. 
Indeed, we will see below that this is not the whole story. 

We have repeated the calculation using the RMF — Dp interaction, i.e. with 
only a p contribution hut tuned to reproduce the same EOS of the RMF — 
{p + 6) case. The results, short-dashed curve of Fig.2, are very similar to the 
ones of the RMF —p interaction. Therefore we can explain the large flow effect 
as mainly due to the different strengths of the vector-isovector field between 
RMF — [p + 5) and RMF — p, Dp in the relativistic dynamics. In fact if a 
source is moving the vector field is enhanced (essentially by the local 7 Lorentz 
factor) [36] relative to the scalar one. 

In order to get the idea we write down, for an idealized situation, the "force" 
acting on a particle. Prom the isovector part of the interaction we get 

5 = ^fS^p^ T f^^p^^ (4) 



for protons and (upper signs), respectively neutrons (lower signs), p* is the 
effective momentum, r the particle proper time (averaged) and pa = Pp—pn the 
isovector baryon density (correspondingly ps3 the scalar one). Here we have 
simplified the problem neglecting the contribution from the current gradient 
in the transverse direction and the current derivative with respect to time. 
We are interested in the difference between the force acting on a neutron 
and on a proton, respectively. Oversimplifying the HIC dynamics we consider 
locally neutrons and protons with the same 7 factor (i.e. with the same speed). 
Then Eq.4 can be expressed approximately by the following transparent form 
(P53 = f^Ps): 



dr dr 



U 
7 



Vp3 



(5) 



where 7 is the Lorentz factor for the collective motion of a given ideal cell. 

Keeping in mind that RMF — (p + 5) has a three times larger p field it is clear 
that dynamically the vector-isovector mean field acting during the HIC is 

much greater than the one of the RMF — p. Dp cases. Then the isospin effect 
is mostly caused by the different Lorentz structure of the "interaction" which 
results in a dynamical breaking of the balance between the p vector and b 
scalar fields, present in nuclear matter at equilibrium. This effect is analogous 
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to the interplay between the isoscalar vector- and scalar-fields which is seen in 
the magnitude and energy dependence of the real part of the optical potential, 
ref. [37]. 




Fig. 2. ^^^Sn+^^^ Sn reaction at 1.5 AGeV (b=6fm): Proton-neutron transverse 
differential collective flow (in MeV/c), vs. rapidity, from the three different models 
for the isovcctor mean fields. Full circles and solid line: RMF — (p + S). Open circles 
and dashed line: RMF — p. Stars and short dashed line : RMF — Dp. 



In order to characterize the effect on differential collective ffows we have cal- 
culated the slope dF'^'^iy) / d{y / yproj) at mid-rapidity. Its value is A&.lMeV/c 
for RMF -{p + 5) and 2'iAMeV/c for RMF - p, i.e. a factor two difference. 

We have also performed some calculations at lower beam energies. We have 
found that up to 500AMeV there is no valuable difference in the differential 
flow predictions among the models discussed here. The effect coming from the 
strength of p field starts to become important around lAGeV, as expected 
from the relativistic mechanism. 

Another interesting observable is the elliptic flow V2{y,Pt), which is derived as 
the second coefficient from a Fourier expansion of the azimuthal distribution 
N{(f),y,pt) — vo{l +ViCos{(l)) + 2v2Cos{2(j))) . It can be expressed as 



>J-z 9 ^ 

Pt 

where pt — ^Jp^. + Py is the transverse momentum [35,38]. 

A negative value of V2 corresponds to the emission of matter perpendicular 
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Fig. 3. Difference between neutron and proton elliptic flow as a function of the 
transverse momentum in the ^^"^Sn +^^2 Sn reaction at 1.5 AGeV b=6fm in the 
rapidity range —0.3 < y/uproj < 0.3. Full circles and solid line: RMF — {p + S). 
Open circles and dashed line: RMF — p. Stars and short dashed line: RMF — Dp. 
The error bars of the RMF — Dp curve are similar to the other cases and are not 
shown just for comprehensibility. 

to the reaction plane, squeeze — out flow. The p^-dependence of f 2, which has 
been recently investigated by various groups [38,39,37], is very sensitive to 
the high density behavior of the EOS since highly energetic particles {pt > 
0.5) originate from the initial compressed and out-of-equilibrium phase of the 
collision, see e.g. Ref . [37] . 

In Fig. 3 we present the pt dependence of the proton-neutron difference of the 
elliptic flow in the same very exotic ^^'^Sn+^^'^ Sn reaction at l.bAGeV (semi- 
central collisions), for mid-rapidity emissions. The larger error bars correspond 
to a reduced statistics when a selection on different pt bins is introduced. The 
effect is increasing for larger pt values due to the smaller number of contribu- 
tions. 

From Fig. 3 we see that, in spite of the statistical errors, in the (p+S) dynamics 
the high-pf neutrons show a much larger squeeze — out. This is fully consistent 
with an early emission (more spectator shadowing) due to the larger repulsive 
p-field. We can expect this appreciable effect since the relativistic enhance- 
ment discussed above is relevant just at the first stage of the coUision. The V2 
observable, which is a good chronometer of the reaction dynamics, appears to 
be particularly sensitive to the Lorentz structure of the effective interaction. 

We have repeated the same set of simulations for the more reahstic ^^'^Sn +^'^^ 
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Sn reaction at 1.5 AGcV (b=6fm), that likely could be studied with the new 
planned radioactive beam facilities at intermediate energies. The results are 
shown in Fig.4. The effect of the different structure of the isovector channel is 
still quite clear, of course with a reduction due to the smaller isospin density 
in the interaction region. Particularly evident is again the splitting in the high 
Pt region of the elliptic flow. 
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Fig. 4. ^'^^S'n + ^^^S'n reaction at 1.5 AGeV (b=6fm) from the three different models 
for the isovector mean fields. Top: as in Fig. 2. Bottom: as in Fig. 3. Full circles and 
solid line: RMF — {p + 5). Open circles and dashed line: RMF — p. Stars and short 
dashed line: RMF — Dp. Error bars: see the text and the previous caption. 

In conclusion intermediate energy heavy-ion collisions with radioactive beams 
can give information on the symmetry energy at high baryon density and 
on its detailed microscopic structure. We have shown that such experiments 
provide a unique tool to investigate the strength of the 5 — like field. The 
sensitivity is enhanced relative to the static property EsymiPs) because of the 
more fundamental covariant nature of the fields involved in HIC dynamics. 

Collective fiows observables are found to be sensitive to isospin effects. Due 
to the time selectivity on the emitted particles the elliptic flow measurements 
appear to be the most appealing, especially for nucleons and light isobars at 
high transverse momentum. 
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